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Abstract 
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1 Introduction 



Currently one of the main interests of nuclear physics is the investigation of hot, 
compressed nuclear matter. One of the most widely used theoretical tools for this 
investigation is the simulation of nuclear collisions by numerical solution of transport 
equations, like e.g. the Vlasov-Uehling-Uhlenbeck (VUU) equation This equation 
and its many relatives have their origin in the kinetic theory of molecular gases. 
Today however they are obtained from relativistic quantum field theory, i.e., by 
processing the Schwinger-Dyson equation for propagators of the interacting quantum 
fields with certain semi-classical approximations |fj, |3|. f4[ . 

It is in principle possible to compare the solution of such transport equations to 
the solution of the full Schwinger-Dyson equation, but to common knowledge such 
an attempt has not been made but once @. The result was a significant difference 
in the time evolution, the full quantum system relaxing somewhat slower than the 
"molecular" transport solution. It has been made plausible that such an effect is 
due to "off-shell" propagation of particles which gives rise to memory effects 

in the hot and dense medium [J7L [|. A derivation of such memory effects within the 
framework of quantum field theory has not been presented so far. 

It is the purpose of the present paper to bridge this gap between the full Schwinger- 
Dyson equation and the traditional "molecular" kinetic picture. 

To this end, one has to go beyond the semi-classical (naive) particle picture when 
processing the integral equation for the propagator. This naive particle picture, tra- 
ditionally labeled "quasi-particle" approximation, assumes the existence of infinitely 
long living elementary excitations even in statistical systems. 

Experimentally, there are indications that this is not very well justified: Spectral 
functions of hadrons in the medium presumably determine the spectrum of emitted 
lepton pairs ||, and "momentum dependent forces" are thought to be important 
for the so-called equation of state in hadronic matter. In relativistic systems such a 
momentum dependence is equivalent to a nonzero imaginary part of the self energy, 
and hence again to a nontrivial hadronic spectral function. Hence one may adopt 
the viewpoint, that the extraction of quantitative information from nuclear collisions 
requires to go beyond the quasi-particle level. 

Also theoretically the naive particle picture is not a good startpoint: Consider the 
breaking of Lorentz invariance in matter or temperature states |10(| . It leads to the 
Narnhofer-Thirring theorem stating the impossibility of a perturbation theory with 
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quasi-particles at finite temperature JTT). In other words, the entities to use for a 
perturbative description in statistical systems must have a continuous mass spectrum 
- they cannot be quasi-particles. How one may formulate a finite temperature per- 
turbation theory for quantum fields with continuous mass spectrum has been shown 
not too long ago [fl^| . Recently we have extended this picture to time dependent 
14] , p~5H as well as spatially inhomogeneous |H| non-equilibrium situations. 

Here, these recent developments will be extended to a general non-equilibrium 
system. The basic assumption for this work is, that nuclear phenomena can be de- 
scribed successfully using methods from relativistic quantum field theory. First, to 
ensure compatibility with existing literature, the Schwinger- Dyson equation for the 
full propagator of an interacting fermionic quantum field is processed according to 
standard procedures. Next, the equation for the spectral function is solved consis- 
tently. The final step of the present paper is a study of the equation resulting from 
insertion of the spectral function into the transport equation. 



2 Matrix valued Schwinger-Dyson equation 

As has been pointed out by various authors, the description of dynamical (time 
dependent) quantum phenomena in a statistical ensemble necessitates a formalism 
with a doubled Hilbert space [0, [|. For our purpose the relevant content of this 
formalism is, that its two-point Green functions are 2x2 matrix- valued. We leave 
it to the reader to chose either the conventional Schwinger-Keldysh, or Closed-Time 
Path (CTP) Green function formalism 0, [13|], or the technically simpler method of 



Thermo Field Dynamics (TFD) (T|, |T| 



Although this choice does not affect the following considerations, it has to be 
stated that the principles of TFD are in accordance with a fundamental require- 
ment: The existence of two mutually commuting representations of the field oper- 
ator algebra. This requirement is not fulfilled by the more conventional CTP (or 
Schwinger-Keldysh) formalism, i.e., CTP may lead to mathematical pitfalls. The 
recently discussed factorization problem seems to point in this direction |B| . 

Within this matrix formulation, we consider the Schwinger-Dyson equation for 
the full fermion propagator of a given model 

S = S + SoQZQS . (1) 
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Here So is the free and S the full two-point Green function of the fermion field, S 
is the full self energy and the generalized product of these is to be understood as a 
matrix product (thermal and spinor indices) and an integration (each of the matrices 
is a function of two space coordinates): 

eS,©g£ = £ /A^Gg. (2) 

3 

Throughout this paper we use the convention to write space-time and momentum 
variables as lower indices, e.g. T, xy = T,(x,y). 

In the CTP formulation as well as in the a = 1 parameterization of TFD [15|, 
the matrix elements of S obey 

S 11 + S 22 = S 12 + S 21 , (3) 

and a similar relation holds for the free propagator Sq. It follows that a linear relation 
exists also among the components of the self energy: 

S ll + y,22 = _ S 12 _ S 21 _ ^ 

Thus, the four components of the Schwinger- Dyson equation are not independent, 
the matrix equation can be simplified by a linear transformation. This linear trans- 
formation, which one may conveniently express as a matrix transformation || 15 1, 
has a physical interpretation only in the TFD formalism: It is the Bogoliubov trans- 



formation defining stable statistical quasi-particle states [16], |17| , not to be confused 
with the traditional quasi-particle approximation. 

In view of this fact, the transformation matrices B will be written in the form 

BW=(f"" ) ""), (5) 

depending on one parameter only. In appendix [A], some useful properties of such 
matrices are listed. Due to the linear relations ((^) and (^) any pair of Bogoliubov 
matrices B(nx), B{n 2 ) gives 

U( \ Q r\ v r\ Q < ft< ^-i (S^Q^ R Q3 R something \ 
%)t 3 5o0S05(%)) =1 S A QE A QS A )' () 

Here, t% = diag(l, —1), S^'" 4 are the retarded and advanced full self energy function, 
and S R,A are the retarded and advanced full propagator (similarly for Sq) 

E H = S 11 + S 12 , S A = S 11 + S 21 

S R = S u - S 12 , S A = S n - S 21 . (7) 
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Together with eqn. (|63| ) follows, that the diagonal elements of the transformed 
equation therefore are retarded and advanced Schwinger-Dyson equation. The off- 
diagonal element is a transport equation. 

We proceed by applying Dirac differential operators to the diagonal parts of the 
simplified equation. They are acting on the free Green function So of a fermion field 
with mass M according to 

S l = \$ x — M => SV S&J = S xy , (8) 

where q = <9^7 M and 8 xy = 5 4 (x — y). The action of these differential operators 
on the off-diagonal elements Sq 2 and Sq 1 of the free propagator gives zero. For the 
retarded and advanced Schwinger-Dyson equation this yields 

l ^xy A = $xy + ^xz A © Sfy A ■ (9) 

To get rid of the dependence on coordinate differences present already in space-time 
homogeneous systems, one performs a Fourier transform of these equations into the 
mixed (or Wigner) representation, according to 

t XP = fd\x - y) exp (iPJx - yY) X xy . (10) 



Here we have set X = (x + y)/2, and the "~" will be dropped from now on. The 
differential operators (inverse free propagators) occurring in the equations ([|) then 
take the form 

^o" 1 = f ~ M + ~$ x =► S^S^ P = 1 • (11) 

The Wigner transform of the convolution S QG as defined in (0) is a nontrivial step: 
Formally it may be expressed as a gradient expansion 

Jd A (x-y) exp(iP M (x-?/r) © G zy = exp (-iO) t XP G XP . (12) 

O is a 2nd order differential operator acting on both functions appearing behind it, 

o = I (d x d G P - d*d°) . (is) 

We have indexed the derivatives to avoid confusion on which factor they act. The 
index E stands for real or imaginary part of the self energy function, and the index 
G stands for action on real and imaginary part of the propagator. Explicitly, the 
first-order term is 

„ I f dG dZ dG \ 
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Obviously, this is the Poisson bracket of the two quantities G and S. The exponential 
of O is a differential operator of infinite order, which formally can be split into real 
and imaginary part as exp (— iO) = cos O — i sin O. 

It is not completely clear how to justify the gradient expansion for arbitrary self 
energy functions, to infinite or to finite order. One may however adopt the following 
viewpoint on the gradient expansion: It is the expansion of the expectation value 
of a composite operator into a series of terms with increasing non-locality Such an 
expansion is always associated with a length scale, and therefore we have to assume 
that the non-equilibrium system is studied only on length and time scales for which 
the gradient expansion is rapidly convergent. 

Since each operator O needs a factor h to render it dimensionless, the gradient 
expansion is also an expansion in powers of h. Truncating it therefore means, that 
quantum effects occurring on very short length or time scales must be absorbed 
into the propagator and the self energy function on larger scales. Indeed, to use 
a fashionable expression for this viewpoint: When taking the gradient expansion 
seriously, we are considering an effective field theory. This requires that the fields we 
are considering have a continuous mass spectrum plf . 

The propagator and the self energy function in mixed representation are split into 
real and imaginary part according to 

Sx'p = G X p =F wA X p 

= S^ + ReS X p T ivrr XP . (15) 

Each of these functions is assumed to be real, and they are spinor-valued. S 5 is 
the part of the self energy function which is local in space and time, which for a 
relativistic model implies that it is the Hartree part of S. 

It is tempting to call A the spectral function of the fermions. However, it cannot 
be guaranteed that the full fermion propagator in non-equilibrium states has a spec- 
tral decomposition [IB]. Therefore, this interpretation has to be treated with caution 
for the moment. 

The gradient expansion and the decomposition of the propagator and self energy 
function are inserted into eqns. (0), leading to 



^-M+i^-exp(-iO)E^^ X p = exp(-iO) 
^-M+i$f x -exp(-iO)£^jG XP = 1 + 



ReSxp Axp + Txp Gxp 



6 



exp(— iO) 



ReZ XP G XP -7i 2 r XP A x fr6) 



The above equations contain a mixture of real and imaginary part, which can be 
separated easily. However, they are not independent in another sense - which can 
be seen when performing the separation and taking the trace over the Dirac indices. 
This gives two pairs of equations: 

Tr [(0 X + 2sinO Y> 5 X ) A XP ] = -2 sin O Tr [Re£ XP A XP + T XP G XP ] 

Tr [($ x + 2smOE s x )G XP ] = -2 sin O Tr [Re£ XP G XP - tt 2 T xp A xp ] (17) 

and 

Tr [(^-M-cosOEy A XP ] = cos O Tr [ReS XP A XP + T XP G XP ] 

Tr [(^-M-cosOEy G XP ] = Tr [1] + cos O Tr [ReT, XP G XP — 7r 2 T XP A X $8) 

The first pair of these equations is obtained from the second pair by application of 
the differential operator 2 tan O. This can be traced back to the fact, that the inverse 
free propagator in eqn. ( |TT1) is a mixture of zeroth and first order in O, formally 



(f-M + ~$ x ) G X p = (l-itanO) (f-M)G XP . (19) 

When the gradient expansion is truncated without considering this fact, this mixture 
may be the source of spurious terms in transport equations [M. 

Obviously, eqn. ( |T8| ) contains only even powers of the differential operator O. 
Solving it in zeroth order O therefore gives a spectral function that is correct up to 
first order in O. 



3 Determination of the spectral function 

We now solve the diagonal components of the Schwinger- Dyson equation, by making 
an explicit ansatz for the Dirac spinor structure of the self energy function: 

Y, x + Re£ XP = S(X, P) + V(X, P) llx + T^(X, P)a, iV , (20) 

where S(X,P), V(X,P) and T(X,P) stand for (Lorentz) scalar, vector and tensor 
component of the real part of the self energy, a^ u = i/2 [7^,7^] and T is antisym- 
metric. To avoid confusion with the symbol used for the fermion propagator, we 
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Figure 1: Two- loop self energy diagrams corresponding to Born scattering diagrams. 



exploit that the scalar and vector piece of the real self energy function only appear in 
combination with M and P. Hence in the following we make use of the abbreviations 



m + s(x, p) p; = p u - v v (x, p) 



(21) 



for effective mass and momentum. Note, that in contrast to other derivations of 
transport equations, the above functions depend on both X and P |20[ . 

Usual derivations of transport equations assume that the self energy contains only 



Lorentz vector and scalar parts f22| , p4| . It is however easy to see that for reasonable 
descriptions of nuclear dynamics a tensor part of the self energy should be taken into 
account. To this end we consider a system, where fermions are minimally coupled 
to bosons. In such a model, the Born diagrams for fermion-fermion scattering are 
obtained by cutting the diagrams in figure p] along the dotted line. If we assume, 
that the Born terms are dominant contributions to the collision term in "molecular" 
transport equations, we thus have to take the diagrams of figure [l] into account in 
our self energy function. 

The contribution of the left-hand diagram in the figure to the fermion self energy 
then contains products of up to three 7-matrices. The threefold product is equiva- 
lent to a pseudovector, which can be excluded in a parity eigenstate. The twofold 
product can be converted to a sum of unit matrix and tensor terms. Thus, if at least 
two independent Lorentz vectors can be specified for a momentum eigenstate, such 



8 



diagrams contribute tensor parts to the self energy functions also for homogeneous 
systems. 

Clearly this is the case for a system at finite density, where the Lorentz invariance 
is broken: Independent Lorentz vectors are the (macroscopic) four-current of 
matter and the momentum of the single-particle state P M (or equivalently the effective 
momentum P*). The tensor part of the self energy is then proportional to the 
commutator [P*^", J^7 M ]. If one wants to consider phenomena like spin diffusion 
20fl , it must not be disregarded. 

By virtue of dispersion relations, a similar decomposition also holds for the imag- 
inary part of the self energy, 



T XP = T S + + r 



t ° iiv 1 



(22) 



It follows from the Dyson equations in the previous section, that the same decompo- 
sition then holds for G and the spectral function: 



Gxp 
A X p 



G s + G^ + Gf <V 



A s + A^ + At a, 



(23) 



The " off-shellness" of the four momentum P u , and correspondingly the generalized 
width of the spectral distribution are given by 

M* 2 



y 

X 

w 



y + n 2 (rj\-r^r£ + 2r fj ^rr) 



2 (P;r£ + M*T S + 2Tp,T?) 
and the propagator denominator is 



V = X' 



7T 



■w 2 



(24) 



(25) 



As was stated after eqn. 
in homogeneous systems 



we may take the solution for A and G as obtained 



It is still correct in first order of the gradient 
expansion, when the local functions for effective mass and momentum are inserted. 
For reasonably smooth systems, or rather on sufficiently large length scales, the 
equation for the spectral function therefore is always an algebraic equation. 



Explicitly, the solutions are [|23 

Gxp = 



f +M* + T^a, u ) X - vr 2 (T s - + If <v) W 



A 



xp 



f +M* + T^a, u W + (r s - + rfV) X 



V 



(26) 
(27) 
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Figure 2: Ghost poles in the perturbative propagator (schematically) 
o are the ghost poles of eqn. (p7|) in a vacuum state, • at nonzero density. ® are the 
particle poles present in certain approximations, like e.g. the one discussed in ref. 
PI). 



These solutions may be obtained in several ways, a straightforward method is sum- 
marized in appendix B. The equations can be combined as 

G XP T ™Axp = {f-M- ReS XP ± mYxp)' 1 
<£> Sxp = (f-M-^fj 1 . (28) 

The inverse on the r.h.s. is to be understood as the proper inverse of a complex 4x4- 
matrix. Hence G and 7rA are real and imaginary part of the same complex function. 
While this may be used to determine the functional form of A in homogeneous states, 
the matrix inversion does not lead to our finding, that A in this form is correct also 
to first order in the gradients. The solution of the differential equations ( [TBI) to this 
order therefore is a necessary step of the derivation (see appendix |B|) . 

However, one may not use the above solutions directly. As is discussed in detail 



in ref. pq| , a perturbative calculation of the self energy function E leads to the 
invalidation of the equation for G. The reason is, that according to Weinberg's 
theorem El the perturbative self energy has a certain asymptotic behavior, which 



leads to the appearance of complex poles in the propagator (|2"8"|), see figure 
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These Landau ghost poles are unphysical, because they appear in the wrong 
complex energy half plane for retarded and advanced propagator - and hence have 
to be removed. In ref. [p8 |, three strategies were discussed for this removal: Direct 
subtraction of the poles, non-perturbative calculation of the self energy function or 
the definition of the full propagator as a dispersion integral 

oo 

G XP t ™A xp = fdE A(X; E,p) ]—— , (29) 

J po - E ± le 

— oo 

with P = (po,p)- It was also shown, that the three strategies are equivalent, and 
that furthermore the non-perturbative calculation of the self energy can be approxi- 
mated by a form factor giving a better asymptotic behavior than the one implied by 
Weinberg's theorem. For the purpose of the present work, eqn. (^) is considered 
the definition of G, in addition to fl26|) . In coordinate space, the above dispersion 
relation implies 

S% A = T2niQ (±(x - y )) A xy . (30) 

For simplicity, consider the spectral function in the case of zero tensor part of the 
self energy. Clearly, this spectral function is not generally of the form 

VV 

with W linear in T. This equation for the fermionic spectral function is only ap- 
proximately correct in case T <C ReE. Such a functional form (which was obviously 
guessed from the non-relativistic result for the spectral function) has been proposed 



A XP = (f + M*\ — 2 (31) 



elsewhere also for arbitrary self energies [^4| , but it is not reproduced by the complete 
derivation carried out here. 

The form (|28| ) of the retarded and advanced propagator implies, how the limit of 
zero T has to be taken, i.e. as 

KmE^ = =Fie 7 ° ) (32) 

to recover the vacuum boundary conditions for the Green functions [p3" |. 

This can be taken into account by adding an infinitely small part to the function 
T. Performing this modification in the above equations then gives for the case of zero 
imaginary part of the self energy: 

f + M* + T^a^ 

Lr) 



*XP 



X 



A XP = (^f + M* + T^a^sign(P*)5(X). (33) 
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In the limit of zero tensor self energy, this turns into 

A XP = [f + M* \ sign(P*) 5 (P* 2 - M* 2 ) . (34) 

Obviously the "particle" -like pole of the Green function is renormalized by the deriva- 
tive of the self energy. 

4 Transport equation 

As was stated above, the off-diagonal component of the Schwinger- Dyson equation 
(in triangular form) is a transport equation. To study this more closely, we specify 
the Bogoliubov parameter of the transformation (||) as 

rii = n 2 = n = const. (35) 

Using eqn. (^) and fl6~4]) , this off-diagonal component then has the form 

S K = S* + Sf H A © S A + S* Y, R © S K - S* E K S A , (36) 

with retarded and advanced components as above and kinetic components 

S? = (1-no) Sf + noS 21 
S K = (1 - n ) S 12 + n S 21 

Ti K = (1 - n ) X 12 + n S 21 . (37) 

As in the preceding section, Dirac differential operators are applied to the off-diagonal 
Schwinger-Dyson equation component. They act on the retarded and advanced prop- 
agator as specified in eqn. (|8|), and annihilate the kinetic component Sq . The kinetic 
component of the matrix valued Schwinger-Dyson equation therefore is 

Together with its adjoint, this equation can be expressed entirely in terms of the off- 
diagonal components of self energy function and propagator. These equations then 



are known as the Kadanoff-Baym equations, see ref. |30 . 

The next step is to transform to the mixed (Wigner) representation, and to insert 
the gradient expansion, 

f - M + ~$ x - exp(-iO) 5&) Si p 

= exp(-iO) [(Re£ XP - ivrr XP ) S$ P + Sf P (G XP + i7rA XP )] (39) 
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The separation of real and imaginary part is a trivial step. However, for comparison 
with existing derivations we perform it explicitly and obtain the two equations 



Tr [ (0 X + 2 sin O (£ x + ReS XP ) + cos O 27rrxp) Sxp\ — 

2iTr [i sin O Sf p G X p - cos O Sf p inA X p] (40) 

Tr [(f - M - cos O (E 5 X + Re£ X p) + sin O nT XP ) S% P ] = 

Tr [cos O Ef P Gxp - i sin O Sf P ivr^p] . (41) 

It is obvious, that these equations are not independent: As was already observed for 
the retarded and advanced equation, ( |4"0| ) follows from (fUD by the application of a 
differential operator. Usually only eqn. fl40f) is processed further, leading to certain 
disadvantages that we will avoid below. 

The final result of the processing will be an equation for the quantity N X p defined 

by 

(1 - N XP ) S% P + N XP fif p = . (42) 
For the purpose of the present paper N XP is taken as a scalar function, it has the 



physical meaning of a generalized phase-space distribution function [31]. The de 



scription of Phenomena like spin diffusion require to use a Dirac matrix valued N XP . 



It follows according to eqns. fl30D and (0), that 



S^ P = 2vri (N XP - n ) A XP . (43) 

Equation ( f42|) has two aspects which are important for the understanding of transport 
phenomena. First, consider it in the equilibrium case: No dependence on X is 
present. The equation is then equivalent to the Kubo-Martin-Schwinger boundary 
condition for the propagator |32| - which implies the proper asymptotic property for 
the quantity N XP we are searching. 

The second aspect of the above equation only makes sense in the context of thermo 
field dynamics (TFD). There, it is the condition for a diagonalization of the full 
propagator by (space-time dependent) Bogoliubov transformations, i.e., according to 
eqn. (§3|) 

B(N XP ) r 3 S XP (S(iVxp))" 1 = ( ° XP ~ mAxP „ t ) ■ ill) 

Gxp + i7T Axp 
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Hence, the full propagator (in mixed representation) is diagonalized by a space- 
time dependent Bogoliubov transformation. One may insert this propagator into 
the Schwinger-Dyson equation - the action of the Dirac differential operator on the 
function N XP would have to be considered then. According to eqn. (|6lD , this would 
only affect the off-diagonal component. The strategy of diagonalizing the full non- 
equilibrium propagator is therefore equivalent to the processing of the Schwinger- 
Dyson equation as done here. 

For the following, we furthermore define a "pseudo-equilibrium" distribution func- 
tion: The property (f|) of the self energy function allows to diagonalize its 2x2 matrix 
structure by a Bogoliubov transformation (see eqn. (|62|) ) with a parameter N xp such 
that 

Y}£ P = 2mN° XP T XP S| P = 2vri (N° XP - l) T XP . (45) 

It is easy to show, that in equilibrium states N° is identical to N as defined in 
([42|) ||15|| . For non-equilibrium states however, N° may vary with time and be quite 
different from N XP . Similar to eqn. ([43]) we obtain 



J XP 



2vri (N° x P - n ) Y 



xp ■ 



(46) 



We now return to the transport equation fl3D|), perform a gradient expansion to first 
order in O and insert the above definitions: 



Tr 



1 - iO) \f - M - + ReS X p) + i7rr XP \ \ {N XP - n ) A XP 



Tr 



1 - iO) | (N° XP - n ) r XP | \g xp + mA XP ^ 



(47) 



The curly brackets contain the quantities acted upon by the O operator according 
to (H). 

The derivative of A on the left side is eliminated using (|17D, and the non- derivative 
terms on the left side are switched over to the right side using (0). All terms 
involving no drop out correctly, and the resulting full quantum transport equation is 



Tr 



A XP (-iO) {f-M-{E 



x 



ReS 



XP 



XP 



N 



XP 



Tr 



-N XP (1 - iO) |r XP | |G X p + wA XP 
h(l-iO) {n xp T xp \ \G XP + mA XP 



(4J 
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5 Vlasov equation 



As the next step, we rederive the "standard" Vlasov equation without collision term. 
To this end, the right side of eqn. ( |4~8| ) is dropped completely. Furthermore, the 
traditional quasi-particle approximation is made: The spectral function is chosen 
according to (0). Note, that according to eqn. (|29"D this is equivalent to an energy- 
independent real self energy function in Wigner representation. In coordinate space, 
such a picture is completely consistent only in an approximation, where all self energy 
functions are local in space and time. In a relativistic model this is the case only for 
the Hartree approximation, nonrelativistically one may also add the exchange (Fock) 
self energy 0. 

In leaving the covariant notation, we set P = (po,P) and X = (t, x). Furthermore 
we assume, that the spectral 5-function (|34"D has at least one pole along the real 
Po-axis, i.e. that one can rewrite it as 

d(P* 2 -M* 2 ) 



S(p**-M* 2 )=J2$ (po ~ Ei(x, P, t) 



dp 



(49) 



where Ei(x,P,t) is the (space-time and momentum dependent) generalized energy 
of the corresponding "particle" -like state. 

Abbreviating the function Nxp at each of the poles as Ni(x,p,t), an integra- 
tion over the energy variable then eliminates the 5-function as well as the energy 
derivative, and one obtains from eqn. (fl8|) 

'dNi(x,p,t) 



E 



at 

dEi (x, p, t) dN t (x, p, t) dEi (x, p, t) dNi (x, p, t) 



. (50) 



dp dx dx dp 

This is the Vlasov (or mean-field) quantum transport equation in standard notation: 
The system evolves in time free of dissipation. 

It is quite instructive to rewrite this in a covariant notation. Using the spectral 
function defined in eqn. fl34|) and y as defined in (p4[), one can transform equation 
([|8|) in a few lines into 

xrv\ l dy dNxp dy dNxp \ n 

This may be integrated over the variable 3^, hence the covariant Vlasov equation 
amounts to the vanishing of the expression in curly brackets on the hypersurface 
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3^ = [[31], pp.438]. The characteristic curves of this partial differential equation are 
parameterized by s, 

ds ~ dPf, ds ~ dX^' ^ ' 

Hence eqn. (pT|) states, that the total derivative of Nxp with respect to s vanishes on 
the hypersurface defined by y = 0. The characteristic equations, together with the 
constraint, are Hamilton's equations for an infinitesimal phase-space "fluid element" 
with proper time proportional to s. In numerical simulations, these two equations 
determine the test-particle trajectories M. 

Note, that the sum over the poles appearing in (|50|) is also inherent to flST]): To 
make the latter meaningful, one has to solve the constraint equation y = 0, and 
in general will find that the hypersurface it defines consists of two or more disjoint 
parts. 



6 Transport equation with quantum effects 



The second step in the study of the quantum transport equation (j48[) consists of 
neglecting only the gradient terms on the right hand side. Real and imaginary part 
separate easily, 

Tr Axp (20) j^p) j !^N XP 

-2 (N XP - N° XP ) Tr [T XP G XP ] ( ->3 ) 
Tr A XP (20) M - (Si + Re£ XP ) j ^N XP 

2tt (N xp - N° XP ) Tr [T XP A XP \ . (54) 

Obviously, this is a transport equation in relaxation time approximation: The deriva- 
tive of N X p is proportional to the difference of N XP and N xp . 

Up to now the spectral function is completely general. However, to study the 
consequences of the off-shell effects in this equation, we have to relate the it as 
closely as possible to the Vlasov equation. This achieved by using the approximate 
spectral function (|3lD, which amounts to an expansion of the complete solution fl27|) 
up to lowest order in T and neglection of the tensor parts. It is inserted in eqn. 
(PJ), i.e., into the imaginary part of the full transport equation. One obtains, after 
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calculation of the trace over the Dirac matrices 

w \ f oy 8N XP ay on xp 



y 2 + vrW J 1 dP^ dX^ dX^ dP» 



2tt W (Nxp - N° XP ) |=0. (55) 



dN XP \ ( dy dN X p dy dN XP 



ds J y=0 \dP, dX» dX, dP» J y=0 v XFJy=0 



While one may consider this equation for each value of po, it obtains a more specific 
meaning when we perform the integration over y. The insertion of the spectral 
function (|31| ) already implies a sufficiently small W, hence the overall factor outside 
the curly brackets can be approximated by a ^-function in this case. 

Then the physical interpretation of W is that of a spectral width multiplied by 
twice the energy of an almost particle-like mode, and the energy of this mode is the 
positive solution of y = (see fl24j)). For this mode we obtain 

2n W (N X P -N XP ). 

(56) 

As argued above s is proportional to the proper time along the characteristic curves of 
the Vlasov equation, and these characteristic curves are the hamiltonian trajectories 
of an infinitesimal phase-space cell, see (l52|) . 

The standard Vlasov equation states, that the occupation probability of this phase 
space cell does not change with proper time. Conversely the equation ( |56|) expresses, 
that such a change of the distribution function occurs in case the spectral function 
peak has a nonzero width. Hence, instead of the Vlasov equation, the result is a 
kinetic equation with dissipation. 

An equation similar to flSUp has been obtained by several authors in spatially 
homogeneous fermionic and bosonic systems |5], [7|, [T5|, [?|, P7fl . For the purpose 

of understanding relaxation phenomena it is of tremendous importance: It implies, 
that without an imaginary part of the self energy function and correspondingly a 
nonzero spectral width there is no relaxation in a quantum system. 

The interpretation of trajectories in phase space is more complicated when pro- 
ceeding to the full spectral function (p7[). Also, when using only eqn. (Q) one obtains 
ugly expressions of the type T s dxM* . . . Hence one has to combine it with the real 
part of eqn. (|4S|), which is written explicitly in (0). The latter then contains terms 
of the type M*d x T s . . . 

The result of this combination is written down only for the sake of completeness: 
It demonstrates how the quantum transport equation including dissipation may be 
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written in covariant form. Using the abbreviations defined in fl2^D , one obtains 
X , , m / 1 



x eqn. (53) + eqn. fl54|) 



]v . ay dN X p _ oy_ on xp \ _ x fdw ox m v o>.v v 



<9P M <9X^ dX^ dP» J \dP^ dX» dX^ dP» 



X 2 f d(X - y) 8N XP d(X - y) dN XP 



tt 2 W V OP,, dX^ dX^ dPv 



2vr (N XP - N° XP ) ((W 2 + X 2 /n 2 ) + 2 (x - (X - y)j } = . 



(57) 



As before, this equation can be studied around a pronounced peak in the spectral 
function. To first order in the peak width, eqn. (|55| ) is recovered. 

Finally we put together all knowledge assembled so far, i.e., the quantum trans- 
port equation ( [48]) is used without further approximation. The crucial problem is 
now to determine the role of the gradient terms neglected in the derivation of fl5T|). 

The right side of fl48|), which contains these terms, has the general structure 



C XP = {1- iO) T XP (G XP + mAxp) ~ exp(-iO) T XP {G XP + mA XP ) . (58) 

r has an obvious meaning to be read off from ([48]) , resembling either T XP or N XP T XP . 

In appendix y of the present paper it is shown, that if and only if the two terms 
on the right side are related through the dispersion relation (|29|) , one may express 
C XP up to first order in the gradient expansion as 

C XP = T XP (G XP + mA XP ) 

+ JdrdEe(-r) exp(-ir£) OT(t + r/2,x;P) A(X; Po + E,p) (59) 

with X = (t, x) and P = (po,P)- Hence the dissipative kernel of eqn. (|5lf ) can be split 
into the time- local piece we have considered already in eqn. (|5"4]), and an integral 
over the past of the system. While the original Schwinger-Dyson equation required 
a time integration reaching into the future of the system, the present equation is 
strictly causal. 

Using the explicit meaning of T from eqn. (^8|), the final result of the present 
work is the equation 



Tr 



Axp (20) If — M — (E s x + ReS X p) } <J N XP 
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2tt (N XP - N° XP ) Tr 



^xp Axp 



+ N XP JdrdE 0(-r) sm(rE) Tr 
- JdrdE e(-r) sin(rE) Tr 



(20) r(t + r/2, f; P) ^(X; p + £, p) 



(20) (iV° m + r/2, x- P) A(X; Po + E, $1)0.) 



The physical interpretation of the integral over the past history of the system is 
obvious: It leads to the memory effects observed in transport theory, and relates 
these to the spectral width of the "particles" propagating in the medium. 

Note, that one may insert a diagrammatic expansion for the self energy func- 
tion, and thereby for T X p and N xp T X p- Consider e.g. the one depicted in figure 
1: In quasi-particle approximation it leads to a standard collision integral on the 
right side of the above equation. However, beyond the quasi-particle approximation 
the oscillating factor sin(rP) may lead to enhancement as well as suppression of 
relaxation. 



7 Conclusions 

In the present work, an equation of motion for the relativistic one-particle phase- 
space distribution function of a fermion system was derived within a real time Green 
function formalism. The only approximation made is an expansion to first order in 
the gradients present in a non-equilibrium system. In particular, no quasi-particle 
approximation was introduced: The full spectral information to the same order in 
the gradient expansion is retained in the dissipative equation fl60D. This places the 
present derivation in the gap between quantum field theory and traditional transport 
equations. 

A result of this study is the observation, that dissipation occurs if and only if 
the spectral function deviates from a 5-function, i.e., if the imaginary part of the 
self energy function is nonzero. For strongly peaked spectral functions, which one 
might associate with "particles" , this implies that dissipation is equivalent to a finite 
lifetime of the excitation. 

Such an interpretation is fully consistent with standard "molecular" transport 
theory: The Uehling-Uhlenbeck collision term in "molecular" transport equations is 
diagrammatically given by Born diagrams ffl. Their closure to a self energy function, 
as depicted in figure 1, gives the lowest order perturbative contribution to the width 
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of the quasi-particle pole in the propagator . This result makes the quasi-particle 
approach to the simulation of nuclear collisions quite dubious. 

Close to equilibrium, the dissipative driving force is proportional to the difference 
between Nxp and N XP , i.e., between the time local distribution function diagonaliz- 
ing the propagator and the one diagonalizing the self energy function. The strategy 
of diagonalizing the full non-equilibrium propagator by thermal Bogoliubov trans- 
formations is therefore equivalent to the solution of the transport equation. One 
may put this into the form of a statement about the TFD formalism: The basic 
idea of thermo field dynamics, the Bogoliubov transformation, achieves the proper 
separation of spectral and statistical information even in a non-equilibrium system 
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To first order in the gradients the spectral function is the solution of an algebraic 
equation, eqn. (P7|). This may be cast together with the definition of the real part 
of the full propagator by dispersion integral, eqn. ( |29| ) and an ansatz to determine 
the self energy from the propagator. Such an ansatz was not discussed in the present 
work, it necessitates a perturbative expansion in terms of the full propagators that 



has been formulated elsewhere [ER [26], E7]. We then have at hand a closed system of 



equations for the Green function of the interacting fermion system, with controlled 
approximations. 

The rewriting of the dissipative kernel in terms of the spectral function and the 
self energy was achieved by using a dispersion integral. This necessitates a proper 
understanding of the analytical structure of the retarded and advanced Green func- 
tion of the system, before following the time evolution of occupation probabilities. 
Hence, the proper ghost-free definition of the real part of the propagator, as outlined 
before and discussed in ref. |2^| , is an essential ingredient of this step. 

The previously made artificial distinction between the spectral width due to col- 
lisions and a "natural" spectral width, like e.g. for resonances [[|, does not arise in 
the present work. However, a distinction is possible between time-local relaxation 
and memory effects present in the system. The influence exerted by the memory 
kernel of eqn. (pCp depends on the details of spectral function and time evolution: 
One may not estimate its effect on the net relaxation time a priori. 

In ref. J7| it was emphasized, that even in a simple model for relativistic heavy-ion 
collisions a slowdown as well as an enhancement may occur through memory effects. 
In H an average slow-down was found in a similar non-relativistic framework, by 
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comparing the solution of the Schwinger-Dyson equation with a "molecular" trans- 
port equation. Accordingly, the next logical step will be a numerical study of the 
derived transport equation floTf), but this is beyond the scope of the present work. 
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A Thermal Bogoliubov matrix 



The thermal Bogoliubov transformations of thermo field dynamics exist in various 



parameterizations. As has been pointed out in ref. [jn| [lf|, the special choice a = 1 
and s = l/21og(l + an) makes them linear in a single parameter n, for fermionic 
fields the explicit form is given in eqn. fl5|). It follows e.g., that 

(ewr^BM = (*(„)>-' ( ° "J) g) B(n) . ( 6 i) 

Obviously, a time- dependence of the Bogoliubov parameter n only affects the off- 
diagonal part of the Schwinger-Dyson equation. 

Now let £ be any 2x2 matrix which fulfills the linear relation (f|). Then one 
obtains 

#(n x ) Sr 3 (B(n 2 )) =1 S A J ' 

Quite similar relations exist for propagator matrices: Let S be any 2x2 matrix 
which fulfills the linear relation (0). Then 

%)r 3 S(%)) =1 I . (63) 

For the purpose of deriving transport equations it is furthermore advantageous to 
exploit a similar relation for the product of three such matrices. Let So also be a 
propagator matrix, whose matrix elements satisfy the relations (|3|). Then one obtains 

(OR yR oR v \ 

S AyA S A J ' ( 64 ) 

where we have used the definition of retarded and advanced components as found in 
eqn. (|?p, and 

x = s^ A s A + s^ R s K 

+ S* ((n 2 - n^S 11 + (n 2 - 1)S 12 - njl 21 ) S A (65) 

with the kinetic components defined in eqn. (p7[). Corresponding relations hold for 
bosonic Bogoliubov transformations, see ref. |L5|]. 



22 



B Fermionic spectral function 



In this section, the equations (|T6| ) are solved using fl20"| ) for the spinor structure of 
the self energy function. First, the trace is taken over the Dirac matrices: 

cos O (P;G£ - M*G S - 2T^Gf ) = 1 - tt 2 cos O (r a A s + T VtlM A% + 2V t ^ u ) 
cos O (P*A% - M*A S - 2T fll/ A? u ) = cosO (T s G s + T v ^G^ + 2T t ^Gn ■ (66) 

Multiplication with 7^ before taking the trace gives the next equations. To eliminate 
the term G^ v one uses fllTf ) , 

cos O (P V *G S - M*G U V ) + 2 sin O (P^Gf - T v » G v>fJ ) 

= cos o (r s A v v + T» v A s ) - 2tt 2 sin O (17 A V4l + T v ^ *4f) 

cos O {P P *A S - M*A V V ) + 2 sin O (P*A^ - T v » A v ^) 

= cos O (r, G u v + V: G s ) + 2 sin O (17 G Vtfl + F Vt(t G?) . (67) 

The remaining spin components are obtained by multiplying with i'j 5 '-f a before taking 
the trace: 

e av ° T cos O (P*G t)(JT - T aT G jV ) = -n 2 e avaT cos O (T t , aT A v ,„ + T v>1/ A t ,ar) 
t oa '° r COS O (P*A t ,ar - T aT Ay^) = e a ™ T cos O {T t , aT G %v + T VjU G t>aT ) , 

where e aucrT is the total antisymmetric tensor in four dimensions. Since X and P are 
independent variables, the sin O-terms from equation fl67D cancel. Furthermore, in 
first order gradient expansion cosO = 1. Reduced spinorial components are intro- 
duced as 

n — S n — V V n — 1 ^ (RSI'S 

AP ~P^~ ~ T 2 ' ^ ' 

similar for g sv t G and k S)V) t <-> T. The scalar equations are 

P* 2 g v -M* 2 g s -2T 2 g t = 1 - vr 2 (AP 2 k s a s + A% + 2T t ^AT) 
P* 2 a v -M* 2 a s -2T 2 a t = M* 2 k s g s + T V>IX G% + 2T t ^ u G? u . (69) 

The vector equations are contracted with P* first: 

9s~ 9v = -7T 2 (k s a v + k v a s ) 

a s — a v = k s g v + k v g s , (70) 
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and then with T, 



p* 2 u n — d v V 
r tv v y s <-J tI i v,v 

P* 2 k a - A V Y 



-Ti 2 {k s A v v T v ^ + a s T v v T v ^) 



9s K v 



v 1 - v,v • 



(71) 



The tensorial equations are contracted with p* u T aT , leading to 

9t - 9v = -k 2 (k t a v + k v a t ) 
a t - a v = k t g v + k v g t , 



(72) 



and with P* y T a t T : 



t,ar l t 

^t,<TT L t 



T 2 g v k t 
T 2 a v k t 



it 2 (a v T t (TT T^ T + k v At.arTI 



9v ^t,ar^t Gt tCrT Y 



(7T 
t 1 



(73) 



We can then use eqns. (^) and (|73|) to eliminate the products of A and G with V. 

I), fE25T), and the final 



The abbreviations used in the following are defined in eqns. 



solution eqn. (|27D is obtained through 



Va s = k s X + W , Va v = -k v X + W , V a t = k t X + 2W 
Vg s = X-k s 7r 2 W , Vg v = X + k v n 2 W , V g t = X - k t n 2 W . (74) 
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C Dissipative kernel of the transport equation 

In this section, the expression on the right side of the quantum transport equation 
( pf ) is evaluated to first order in the gradient expansion. First, we write it to infinite 
order in O 

C XP = exp(-iO) f XP {G X p + inAxp) ■ (75) 

The next step consists of using (|30|) in connection with ( |12|) (See remark in the text 
after eqn. (|5Tf)). We remove the gradient expansion for this intermediate step, 

C XP = 27rJd 4 z^e(-z )e W (^(P-Q)z+ 1 -d x dij x 

f (X + z/2- P) iA{X; Q) . (76) 

The three-dimensional integrations can be carried out, with X = (t,x) and P = 
(po,p) one obtains 

C X p = i JdrdE 6(-r) exp ^-iO - \rE + i $d£j X 

f (t + r/2, x ; P) A(X; Po + E,p). (77) 

The solutions for G and A are known only to first order in the gradient expansion, 
hence this expression is also valid only in first order in O. Furthermore, the additional 
derivative term in the exponent is eliminated by shifting the argument of V: 

C X p = i JdrdE Q(-t) exp(-irE) (l -iO+ ^d£\ x 

f(t + r/2, x ; P) A(X; p + E,p) (78) 
= i JdrdE 0(-r) exp(-ir£) (f (X ; P) A(X ; p + E , p) 

-iOf(t + r/2, x; P) p + E,p) 

In the first part, the i- integration can be carried out, 

Cxp = fdE 1 ——f(X;P)A(X;E,p 

J po-E-ie 

+ JdrdEQ(-r) exp(-irE) of (t + r/2, x; P) A(X;p + E,p) . (79) 
Together with (|29|) , this is equivalent to eqn. (BH). 
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